Why does wurtzite form in nanowires of III-V zinc-blende semiconductors? 
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We develop a nucleation-based model to explain the formation of the wurtzite (WZ) phase during 
the vapor-liquid-solid growth of free-standing nanowires of zinc-blende (ZB) semiconductors. Nucle- 
ation occurs preferentially at the edge of the solid/liquid interface, which entails major differences 
between ZB and WZ nuclei. Depending on the pertinent interface energies, WZ nucleation is favored 
at high liquid supersaturation. This explains our systematic observation of ZB during early growth. 

PACS numbers: 68.65.La,64.60.Qb,81.05.Ea,81.15.Kk,64.70.Nd 



Free-standing wires with diameters ranging from hun- 
dreds down to a few nanometers are nowadays commonly 
fabricated from a large range of semiconductor materials 
[E S 0, IH HI • These nanowires (NWs) have remarkable 
physical properties and many potential applications. The 
present work deals with the epitaxial growth of NWs of 
III-V semiconductors on a hot substrate. Metal catalyst 
nanoparticles deposited on the substrate before growth 
define the wire diameter. According to the vapor-liquid- 
solid (VLS) growth mechanism, the atoms are fed from 
the vapor phase to the solid wire through this particle 
(or droplet), which remains liquid during growth Q. 

We consider III-V compounds which, under bulk form, 
adopt the cubic zinc-blende (ZB) crystal structure 0] (al- 
though some non-ZB high-pressure phases Q may be 
metastable at atmospheric pressure 0]), leaving aside 
nitrogen-based NWs. We discuss the usual case of NWs 
grown on a [111]B (As-terminated) face of the ZB sub- 
strate. Probably the most surprising feature of these 
NWs is that, in contrast to their bulk counterparts, they 
often adopt the hexagonal wurtzite (WZ) structure. This 
was observed for most ZB III-V materials and growth 
techniques [USB E) EH- However, although often 
dominantly of WZ structure, the NWs usually contain 
stacking faults (SFs) and sequences of ZB structure. The 
coexistence of two phases is clearly a problem for basic 
studies as well as applications, so that phase purity con- 
trol is one of the main challenges of III-V NW fabrication. 

The surprising prevalence of the WZ structure in III-V 
NWs has not been explained satisfactorily so far. Here, 
based on new experimental observations, we propose 
an explanation of the occurrence of the WZ structure 
and develop a model predicting quantitatively in which 
growth conditions it should form. We consider the spe- 
cific case of gold-catalyzed GaAs NWs grown by molec- 
ular beam epitaxy (MBE) on a GaAs substrate but we 
expect our model and our conclusions to remain valid for 
any ZB III-V compound and any growth method. 

Let us start with briefly reviewing previously proposed 
explanations. Calculations give the difference Sw in co- 
hesive energy between ZB and WZ bulk GaAs as about 
24 meV per III-V pair at zero pressure Q. It has been 
argued that this favoring of the ZB form might be offset 



in NWs of small diameter by the large relative contribu- 
tion to the total energy of either the lateral facets [T3| 
or the vertical edges separating the latter [3] (provided 
the specific energies of these features are less for WZ than 
for ZB). This naturally leads to critical radii under which 
WZ NWs should be more stable than ZB NWs. For in- 
stance, Akiyama et al. calculated a critical radius of 5.6 
nm for GaAs NWs. These approaches have in common 
to treat the energetics of fully formed NWs and to pre- 
dict critical radii far too small to explain the occurrence 
of WZ in NWs with radii up to at least 100 nm. 

On the other hand, from the very beginnings of VLS 
studies, it has been argued that the two-dimensional (2D) 
nucleation of new solid layers from the supersaturated liq- 
uid was of paramount importance [3] and most theories 
of NW growth take nucleation into account 15, 3, TtJ • 
The fact that the faults in each phase and those separat- 
ing ZB and WZ regions are perpendicular to the growth 
axis, in other words that each monolayer (ML) of III- 
V pairs is uniform in structure and orientation, strongly 
suggests that, once a nucleus of critical size is formed, it 
rapidly spreads out laterally over the whole solid/liquid 
(SL) interface [B|, [ll|, unless the wire is very wide . If 
so, the reason for the formation of the WZ phase should 
not be searched in the energetics of fully formed NWs 
lil lH but instead in the preferential formation of WZ 
nuclei at the SL interface. This is the aim of the present 
work. 

Given the prominence of WZ in GaAs NWs, it is in- 
teresting to find two instances where the ZB structure 
systematically appears. The first one corresponds to the 
initial stage of NW growth. By growing GaAs NWs for 
short times, we ensured that the foot of each NW could 
be observed clearly before being buried by the 2D layer 
which grows concomitantly albeit more slowly between 
the NWs. The cross-sectional images obtained ex situ 
by transmission electron microscopy (TEM) prove that 
the whole NW is initially pure ZB. Growth then switches 
abruptly to WZ stacking (Fig. [T]) . Scanning electron mi- 
croscopy shows that, at this early stage, the NWs are 
pyramids with triangular bases and tilted lateral facets. 
These become vertical at the ZB^WZ transition (Fig.[l|. 

The second observation is that when we terminate 
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FIG. 1: (Color online) TEM image of a short GaAs nanowire 
with high resolution close-up of the ZB— >WZ transition zone. 



MBE growth by switching off the Ga flux while maintain- 
ing an As flux, a section of NW grows that systematically 
adopts the ZB structure. This effect, already reported for 
GaAs NWs grown by chemical beam epitaxy [3J] , is inter- 
preted as a partial consumption of the Ga dissolved in 
the gold particle to form the terminal section of the NW. 

These two situations where ZB forms have in common 
to be transient growth phases during which the supersat- 
uration of Ga (and possibly As) in the liquid is less than 
during steady NW growth. Before growth, the deposited 
Au droplets dissolve the substrate locally to achieve equi- 
librium with it: the supersaturation is zero. When vapor 
fluxes are turned on, the supersaturation increases until a 
permanent regime settles. Symmetrically, during growth 
termination, the Ga concentration in the droplet, and 
hence supersaturation, decrease, since the atoms used to 
build the NW are not replaced. This strongly suggests 
that ZB systematically forms when the supersaturation 
is less than some critical value and, conversely, that WZ 
formation requires a high supersaturation. 

This confirms the importance of nucleation. Indeed, 
according to nucleation theory, the work needed to form 
solid nuclei from a fluid phase is maximum for a critical 
nucleus size flj| . If this energy barrier is less for a certain 
crystal structure than for another, the first may nucleate 
preferentially even if the second one is more stable in 
bulk form. Since critical sizes and energy barriers scale 
inversely with supersaturation, our observations point to 
a preferential formation of WZ when the critical nuclei 
are small. This recalls the preference for WZ of NWs of 
small radius 
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hold for nuclei. Instead, we should compare the probabil- 
ities of forming various nuclei from the liquid phase. As 
a first approximation, we shall compare their formation 
energies, which largely determine these probabilities. In 
order not to obscure our demonstration, we keep as far 
as possible to continuous nucleation theory (CNT) [l8j |. 
resorting to an atomistic picture only when necessary. 
We proceed in two stages. We first show that nucleation 
must occur at the edge of the SL interface (the triple 
solid/liquid/ vapor line) rather than elsewhere in the SL 
interface. We then show that, along this line, the forma- 
tion of WZ nuclei may be favored over that of ZB nuclei. 





FIG. 2: (Color online) (a) Nucleus at the NW/liquid interface, 
with interfaces of interest, (b) Transferring the nucleus to the 
triple line eliminates and creates interfacial areas (thick lines) . 



Let us consider the interface between a NW (hereafter 
called substrate) and the liquid droplet. According to 
CNT, forming there a solid epitaxial 2D island of height h 
(I ML), perimeter length P and upper area A (Fig.[^a)) 
from the liquid phase, involves a change of free enthalpy: 

AG = -AhAn + PhjiL + A (j NL - j SL + Isn) (1) 

In Eq. (JlJ, A/i > is the difference of chemical potential 
for III-V pairs between liquid and solid phases, per unit 
volume of nucleus; jil is the energy per unit area of the 
lateral interface between nucleus and liquid; 7sl, 7s at 
and 7tvl are, respectively, the energies per unit area of 
the substrate/liquid, substrate/nucleus (SN) and upper 
nucleus/liquid (NL) interfaces (Fig. (Ha)). 

A given nucleus (set of atoms with fixed relative posi- 
tions) of ML height cannot be said to be of ZB or WZ 
structure. It is only the orientational positioning of the 
nucleus with respect to the previous ML which deter- 
mines if the stack of 2 MLs formed by adding the nucleus 
is of the type found in ZB or WZ crystals (Fig. [3]). In the 
former case (hereafter 'ZB position') the GaAs4 tetrahe- 
dra have the same orientation if the Ga atom belongs 
either to nucleus or to previous ML whereas tetrahedra 
and nucleus are rotated by an odd multiple of 7r/3 in the 
latter case ('WZ position') Q. ZB and WZ sequences re- 
quire the nucleation of each ML in, respectively, ZB and 
WZ position with respect to the previous ML. 

Whatever the position of the nucleus, interfaces SL and 
NL have the same atomic configuration (Fig. [3]) so that 
Isl = Inl- On the other hand, nucleation in WZ posi- 
tion costs some SN interfacial energy ('jYn — If > 0) 
whereas nucleation in ZB position does not (j§§ = 0). 
Since forming a WZ sequence on a ZB substrate creates 
a SF, is sometimes taken [j| as half the SF energy in 
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FIG. 3: (Color online) A given nucleus (open symbols) in ZB 
(left) and WZ (right) positions on top of a (lll)B ML (disks). 
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the bulk ZB phase [l9| ■ Finally, from Eq. (P) , the forma- 
tion enthalpies of a given nucleus in ZB or WZ position 
are AG ZB > WZ = —AhAfi + Ph ltL + A-/^ wz . Since 
AG WZ - AG ZB = Ay F > 0, ZB nucleation is favored. 

In the foregoing discussion, we did not specify if jil 
refers to ZB or WZ because the nucleus, which clearly has 
the same edges in both positions, was assumed to be lat- 
erally surrounded by the liquid. However, we now show 
that nucleation should take place at the triple line. Let 
us compare the formation of a given nucleus at two differ- 
ent locations (Fig. [2]): its lateral surface is either entirely 
surrounded by the liquid (1) or partly surrounded by the 
vapor because of nucleation at the triple line (2). We now 
have to distinguish the specific energies ^il and ytv of 
the lateral nucleus/liquid and nucleus/vapor interfaces. 
The key point is that shifting the nucleus from location 
(1) to location (2) at constant liquid volume has a major 
effect: it eliminates part of the liquid/vapor interface and 
replaces it by nucleus/vapor interface (Fig. [2]). Suppose 
that forming an area s of nucleus/vapor interface elimi- 
nates an area rs of liquid/vapor interface, of energy jlv- 
If a is the fraction of the island perimeter in contact with 
the vapor, the formation enthalpy of the nucleus now is: 

AG = -AhA/j, + Ph[(l-a)-f iL + a(j ev -t^/lv)] 
+Aj SN (2) 

where, as above, jsn — or 7^ for, respectively, ZB 
or WZ nuclei. For a given nucleus, the difference in 
formation enthalpies between locations (1) and (2) is 
AG2 — AGi = aPh (7fy — 7^l — tjlv)- Factor r cannot 
be calculated exactly. We estimate it by considering an 
artificially axisymmetric nucleus with a vertical lateral 
surface contacting the vapor along the whole triple line. 
This geometry preserves a spherical liquid/ vapor inter- 
face. It readily yields r = sin j3, with (3 the contact angle 
between droplet and substrate (Fig. [2|) . In all our sam- 
ples, 90° < < 130° (after growth) so that 0.85 < r < 1. 

Hence, a given nucleus tends to form at the triple line 
if Jev — Iil — tjlv < 0, with t ~ 1. Before growth, the 
contact angles between our liquid droplets and bulk GaAs 
are close to 7r/2, which implies (from Young's equation) 
that the solid/liquid and solid/vapor interface energies 
are close for (lll)B surfaces. Assuming that this also 
holds for the lateral nucleus faces yields ju, ~ Jev- I n 
turn, 7lv should lie between the surface energies of pure 
liquid Au and Ga (1.14 and 0.72 J.m" 2 0). Hence, the 
above inequality is safely satisfied and the critical nuclei 
should form at the edge of the droplet. In short, form- 
ing the nucleus there is advantageous because it elimi- 
nates a portion of the pre-existing droplet surface; this 
largely outweighs the replacement of part of the lateral 
nucleus/liquid interface by a possibly slightly costlier nu- 
cleus/vapor interface. Note that nucleation at the triple 
line in GaP NWs has previously been argued for on an 
entirely different basis, namely the low solubility of phos- 



phorus in gold [5j. Our argument is of general validity 
and would even apply to solid catalyst particles 

To demonstrate the advantage of WZ over ZB nucle- 
ation at the triple line, we consider the atomic structure 
of the edges of the top facet of the NW and of the nucleus. 
We simply select the low energy configurations discussed 
below by restricting ourselves to at most one non- vertical 
dangling bond per atom and to stoichiometric nuclei. At 
the first ZB— >WZ transition (Fig. Q]) , given the pyramidal 
shape of the NW foot, the three edges of the top facet 
must be rows of Ga atoms. A given nucleus can again 
be placed there in ZB or WZ position, but this now en- 
tails major differences. The ZB nucleus must itself have 
a Ga edge at the triple line, so that the edges of the top 
NW facet and of the nucleus develop a tilted {111}A- 
type lateral facet (Fig. HJa)). The WZ nucleus has an As 
edge at the triple line; the lateral facet is then vertical 
and of {1010} WZ type (Fig.^b)). Since WZ sequences 
require the repeated nucleation of each ML in WZ posi- 
tion, the latter must therefore be favored at both Ga and 
As edges. This is also needed when WZ NWs with six 
{1010}-type vertical facets have started developing since 
their top facets have three edges of each type. However, 
at an As edge, ZB and WZ positions also differ: in ZB 
position, the nucleus itself has an As edge and generates 
a {lll}B-type tilted lateral facet (Fig. HJc)), whereas in 
WZ position (not shown) it has a Ga edge and the lateral 
facet is again vertical, of {10l0} WZ type. 

Let us first compute the energy changes Aj upon trans- 
ferring a given nucleus from location (1) (Fig. Ela)) to 
three possible locations at the triple line (Fig. [4]), one in 
WZ position with a lateral {1010} facet (j — W) and two 
in ZB position with lateral {111}A (j = A) or {111}B 
(j = B) facets. When the lateral facets are tilted toward 
(j = A) or away from (j = B) the NW axis, the transfers 
also modify the areas of the solid/liquid and liquid/vapor 
interfaces, at fixed liquid volume. This can again be 
estimated by first considering toroidal nuclei preserving 
spherical liquid/ vapor and axisymmetric SL interfaces, 
and then taking Aj as proportional to the fraction aP of 




FIG. 4: (Color online) Top (top row) and side (bottom row) 
views of a given nucleus in ZB (a) and WZ (b) positions at a 
Ga edge and in ZB position at an As edge (c). Edges are in- 
dicated by arrows (top), the vapor side by V and non- vertical 
dangling bonds by segments (bottom). Key as in Fig. [3] 
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the perimeter of the actual nucleus in contact with the 
vapor. We find Aj = aPh (7,- — jil — Jlv sin [3) where 
7j = Jj/ cos Oj + (j LS + "f LV cos (3) sinOj, with 7, the in- 
terface energy between vapor and lateral facet and 9j the 
angle of the latter with the vertical (6w — 0,6b — ~0a — 
19.5°). Our 7j have the same expression as the energies 
calculated by Ross et al. in a 2D model of facetted NWs 
2 if . but they pertain to ML high nuclei. 

Although it is not essential, we now specify that the 
nuclei are equilateral triangles of side D, one of them at 
the triple line (a — 1/3). Their formation enthalpies are: 

/3 

AG j = -^-D 2 hAn + ZDh liL -Dh{ liL + lLV sinl3) 



2 



V3^ 2 



D hj SN 



(3) 



must be high enough, in agreement with our experimen- 
tal results. Our aim was to identify important effects and 
parameters, not yet to give a complete description of the 
complex interplay of the two phases. We now intend 
to calculate the actual nucleation probabilities (includ- 
ing the effects of temperature and geometry), evaluate 
more precisely the energies of various nuclei (including 
non-stoichiometric ones) forming on NWs with different 
cross-sections, and take into account growth conditions 
in more details, in particular the supersaturation of each 
atomic species which appears here only indirectly via sur- 
face energies. 

This work was partly supported by the SANDIE Net- 
work of Excellence of the European Commission (Con- 
tract No. NMP4-CT-2004-500101). 



Only the last two terms differ between nuclei and the 
last one is non-zero only in WZ position. Maximizing 
AGj with respect to D yields the critical size D* and 
the energy barriers AG* = AGj (£>*) for each nucleus: 



AG* W = 



V3 



6 Afj, 
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AGt = — -± for k = A,BU) 
k 6 A^i ' 



where T, 



j — 2j£L + 7j — 7Lysin/3 is an effective edge 
energy. WZ nuclei dominate if AG*y < AG*, for k = 
A, B. This requires two conditions. The first one, 7n/ < 
jk for k = A, B, is material-related and involves only 
interface energies. A second, growth-related, condition is 
that the supersaturation be larger than a critical value, 



Afi* = max fe= ^ s 



IF \ 



to overcome the SF. 



As a first approximation, we estimate the ZB 7^ ener- 
gies from those of (111) A and B surfaces computed for 
As-rich (MBE) vapors, namely 0.82 and 0.69 J.m" 2 f^. 
In the extreme cases of drop surfaces of pure Ga (sur- 
face segregation of the low energy atom) and pure Au, 
this gives respectively 74 = 0.76 and 0.83 J.m -2 and 
7b = 0.84 and 0.77 J.m~ 2 , well above the low 7111s en- 
ergy. Conversely, 7^ = 71010 is unknown for As-rich va- 
pors. According to the previous discussion, WZ forms be- 
cause 71V < 74 and 73 ■ This hypothesis is strengthened 
by calculating the critical supersaturations for a plausi- 
ble range of such values, 0.7 < 7V < 0.75 J.m -2 . For 
e.g. a Ga drop surface and (3 = 120°, A/i* ranges be- 
tween 230 and 1570 meV, which is indeed of the order of 
our experimental supersaturations (several 100 meV). 

In summary, we developed a nucleation-based model to 
explain the occurrence of the WZ phase in nanowires of 
ZB semiconductors, at least at certain stages of growth. 
A key and general result is that 2D nucleation takes 
place preferentially at the edge of the solid/liquid inter- 
face. When formed at this triple line, WZ and ZB nuclei 
present major differences and WZ nucleation is actually 
favored for certain ranges of the interface energies in- 
volved. In addition, the supersaturation of the liquid 
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